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ABSTRACT 

We present empirical scaling relations for the significance of absorption fea- 
tures detected in medium resolution, far-UV spectra obtained with the Cosmic 
Origins Spectrograph (COS). These relations properly account for both the ex- 
tended wings of the COS line spread function and the non-Poissonian noise prop- 
erties of the data, which we characterize for the first time, and predict limiting 
equivalent widths that deviate from the empirical behavior by < 5% when the 
wavelength and Doppler parameter are in the ranges A = 1150-1750 A and 
b > 10 kms -1 . We have tested a number of coaddition algorithms and find the 
noise properties of individual exposures to be closer to the Poissonian ideal than 
coadded data in all cases. For unresolved absorption lines, limiting equivalent 
widths for coadded data are 6% larger than limiting equivalent widths derived 
from individual exposures with the same signal-to-noise. This ratio scales with 
6-value for resolved absorption lines, with coadded data having a limiting equiv- 
alent width that is 25% larger than individual exposures when b f» 150 kms -1 . 

Subject headings: line: profiles — methods: data analysis — ultraviolet: general 



1. Introduction 

The Cosmic Origins Spectrograph (COS) was installed aboard the Hubble Space Tele- 
scope (HST) in May 2009 as part of HST Servicing Mission 4. COS is a slitless spectrograph 
with far-UV (FUV) and near-UV spectr oscopic and imagin g; mo des. Its design philoso phy 



and on-orbit performance are detailed in iGreen et al.l ( 120121 ) and lOsterman et al.l (120111 ) 



Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Associated Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 
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COS has three FUV gratings: the medium resolution G130M and G160M gratings, 
which use five central wavelength settings to cover the wavelength ranges 1135-1470 A and 
1385-1795 A, respectively, at a spectroscopic resolving power of R = A/AA w 18000; and the 
low resolution G140L grating, which uses two central wavele ngth settings to cover the wave- 
lengt h range 900-2150 A at a resolving power of R ~ 2500 (IDixon et al.ll2011t iGreen et al. 



20121 ). Three additional c entral wavelengths for the G130M grating have become available re- 
cently ( IDixon et al.ll201lh . but since they are not as well-characterized as the original G130M 
settings we do not investigate them here. All of the gratings disperse light, re-image the di- 
verging telescope beam, and correct for the astigmatism of the spectrograp h design and the 



spherical aberration of t he HST primary mirror in a single optical element (lOsterman et al. 



201 ll ; IGreen et al.l 120121 ) . This design minimizes reflective losses in the system, creating a 
spectrograph with unprecedented FUV sensitivity; the COS FUV thro ughput is 10-20 tim es 
greater than comparable modes of previous HST FUV spectrographs (IGreen et al.ll2012l ). 



Thermal vacuum testing of the COS resolving power pre-launch found that the line 
spread function (LSF) of the G130M and G160M gratings was well-approximated by a Gaus- 
sian with a FWHM of ~6.5 pixels (At> ~ 17 kms -1 ). However, the measured on-orbit LSF 
of these (and all other) gratings was found to have extended non-Gaussian wings that vary 
as a function of waveleng th. These wings are cau sed by mid-frequency wavefront errors in 
the HST primary mirror (jGhavamian et aL 2009 ) and are further enhanced by a small but 
measurable amount of scattering caused by microroughness on the primary mirror surface 
(jKrissll201ll ). These deviations from a Gaussian profile affect COS more than HSTs previous 
UV spectrographs because of its slitless design, and reduce the instrument's sensitivity to 
weak absorption features by decreasing the resolution and intensity at the core of the LSF. 



In this Paper, we expand upon the work of iGhavamian et al.l (120091 ) and iKrissI (1201 ll ) by 
generalizing their equivalent width formalism to spectroscopically resolved lines and investi- 
gating the effects of non-Poissonian noise in COS G130M and G160M data on the significance 
of absorption features. While important for all scientific applications, this non-Poissonian 
noise affects the detectability of broad, shallow absorption features in particular , such as Lyq 



and metal- l ine absorption associated with the warm-hot intergalactic medium (|Savage et al. 



201(1 1201 it iNaravanan et al.ll201ll : iDanforth. Stocke. fc Shullll2010bl : iDanforth et al.ll201ll ). 
In Section [2] we demonstrate the presence of non-Poissonian noise in COS FUV data and 
characterize its effects on the limiting equivalent widths of absorption lines. In Section [3] we 
derive empirical approximations for the significance of absorption features detected in COS 
FUV data that properly account for both the extended wings of the on-orbit COS LSF and 
the non-Poissonian noise properties of the data. These approximations predict equivalent 
widths that differ from the empirical behavior by < 5% over the entire G130M + G160M 
bandpass (1150-1750 A) when the Doppler parameter b > 10 kms -1 . Finally, we summarize 
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our results and most important conclusions in Section HI 



2. Non-Poissonian Noise Properties 



When examining a spectrum, one is often interested in knowing the smallest equivalent 
width that can b e measured to a c e rtain significance level. This limiting equivalent width 
was expressed by iGhavamian et al.l (120091 ) as 



N a \ x 1 



lim 



(S/N). f e (x) 



X 



where x is the width (in pixels) of the discrete region of integration over which the limit is 
calculated, N a is the significance level of the limit expressed as a multiple of a for Gaus- 
sian distributed errors, X x = AX x x is the width of the integration region in wavelength 
space, AA is the dispersion o f the spectrum (9 .97 mA pixel" 1 for the G130M grating and 
12.23 mA pixel" 1 for G160M; Dixon et al.ll201lh . (S/N) a is the signal-to -noise ratio of the 
spectrum at a resolution of X x , and f c (x) is the fractional area of the LSF contained within 
the region of integration. 

If the noise in the spectr um is purely Poissonian, (S/N) x = (S/N) 1 xx 1//2 , where (S/N) 1 is 
the signal-to-noise per pixel. IGhavamian et ajj (120091 ) use this relation to derive the limiting 
equivalent width to be 

N a AX x 1/2 
(S/N) 1 f c (x) 

The limiting equivalent width is obviously a function of x, the size of the region of integration. 
Fortunately, x 1 ^ 2 / f c (x) can be uniquely minimized such that there is an optimum integration 
width, x opt , that provides the most stringent equivalent width limit. Note that for COS 
data f r is a function o f both x an d A since the wings of the LSF vary with wavelength 
( IGhavamian et al.l 120091 ; iKrisd 1201 ll ). Thus, x opt and WW are wavelength dependent even 
for fixed signal-to-noise. 

However, the noise in COS FUV data is not purely Poissonian. We define the measured 
relationship between the smoothed signal-to-noise and the signal-to-noise per pixel to be 



T](x) = 



(S/N), 
(S/N), 



(3) 



where T](x) = x 1 ^ 2 for purely Poissonian noise. To quantify the functional form of rj(x) 
we have measured (S/N), as a function of x at several wavelengths for all of the FUV 
data obtained as part of COS GTO programs 11520 and 12025 (G130M + G160M spectra 
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of 16 QSOs with coadded (S/N) 1 ~ 5-25). All exposures were reduced with CalCOS 
v2. 17.3a and proc essed with cust o m IDL rout ines for aligni n g and coadding COS FUV 



data as detailed in iDanforth et al.l (l2010al ) and iKeeney et al.l (J2012J). The signal-to- noise 
ratio was calculated by identifying line-free continuum regions, smoothing the spectra by x 
pixels, and measuring the rms continuum deviations in the smoothed spectra. 

Figure [1] shows our measured values of rj(x) as a function of x for both individual 
exposures (diamonds) and coadded data (asterisks). The plotted symbols represent the 
average value of 77 as determined from our 16 sight lines, each of which had measurements 
taken at 13 different wavelengths. We find that 77 is well fit by a power law of the form 

T](x) = 7] + X 13 . (4) 

The solid line in Figure [TJ indicates our best-fit power law for coadded data, and the dot- 
dashed line indicates our best-fit power law for individual exposures. The dashed line il- 
lustrates the behavior of Poissonian noise. Only smoothing lengths larger than the FWHM 
of the on-orbit LSF (dotted vertical line; see Table [2]) contribute to the fits. Our best-fit 
values of r] and (3 for Poissonian data, individual exposures, and coadded data are listed in 
TableU 

We find a significant difference in rj(x) measured from individual exposures and coadded 
data, but emphasize that even though the non- Poissonian noise properties of coadded data 
are more severe than for individual exposures, the coadded signal-to-noise is significantly 
higher, allowing observers to set more stringent equivalent width limits from coadded data 
in most cases (see discussion in Section \3.1\i . We find no evidence for systematic trends in 
rj(x) as a function of wavelength or signal-to- noise in either individual exposures or coadded 
data. For coadded data, we also find no evidence for systematic trends in T](x) as a function 
of the number of exposures being added together. 

The G140L grating is not the best choice for dete cting weak spectral f eatures, so we 
do not describe it in detail here. However, we note that ISyphers et al.l (120121 ) find rj ~ x 0A4 



for small smoothing lengths (x < 14 pixels), and 77 ~ x 0A2 for larger smoothing lengths 
(x < 40 pixels) in individual G140L exposures. These values are both close to our measured 
value of rj ~ 0.08+x ' 42 for individual G130M and G160M exposures (see Table[T]), suggesting 
that the non- Poissonian noise in individual COS exposures originates in the FUV detector 
itself rather than any particular grating. 

The COS FUV detector utilizes two microchannel plate (MCP) segments, which both 
have dead spots, gain variations, and shadows from the ion-repeller grid wires used to improve 
detector quantum efficiency. Additionally, each segment has irregularities introduced by 
hexagonal and moire patterns in the MCPs themselves. All of these features create fixed- 
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pattern (i.e., non-Poissonian) noise in COS FUV spectra when these spatially variable two- 
dimensional "pixels" are compressed to one dimension as part of the spectral extraction. 
While some of these effects (e. g., grid wire sha dows) can be corrected by flat fielding, they 
cannot be removed altogether (jAke et al.ll2010l ). 



Standard observing strategy suggests that targets should be observed using multiple 
FP-POS positions and / or central wavelengt h settings to reduce the effects of fixed-pattern 
noise (see Section 5.8.2 of iDixon et al.ll201ll ). but we find just the opposite effect. Figured] 
shows that r](x) for coadded data is further from the Poissonian expectation than r)(x) for 
individual exposures. This result is not in itself surprising given the origin of the non- 
Poissonian noise in individual COS exposures and the spectrograph's slitless design, since 
the hexagonal and moire patterns of the MCPs can create large fixed-pattern differences 
from small variations in the spectrum's position in the detector's cross-dispersion direction. 
In general, these differences do not simply average out, and can potentially constructively 
interfere to create amplified fixed-patte rn noise; indeed , this is one of the reasons why a 
robust one- dimensional FUV flat field (jAke et al.l 120101 ) has been difficult to produce for 
COS. 

If the signal-to-noise of a dataset is high enough (S/N > 10 per pixel in individual expo- 
sures) an iterative process can be utilized, which was first developed for high signal-to-noise 
observations obtained wit h the Goddard High Resolution Spectrograph (ICardelli &: Ebbets 
19931 ; lLambert et al.lll994j ). This iterative procedure can be used to create a custom flat field 
that improve s the final signal- to-noise of individual exposures and coadded data for that spe- 
cific dataset ( jAke et al.ll2010l ). Given the higher signal-to-noise achieved in data processed 
by these custom flat fields compared to the standard flat fields applied by CalCOS, the 
non-Poissonian behavior shoul d be less severe when a custom flat field is employed. The 
Mrk 421 sight line studied by iDanforth et al.l (120111 ) is the only one in our sample with 
sufficient signal-to-noise to utilize this technique, and we note a modest improvement in 
the noise properties of this dataset when a custom flat field is employed (as was done in 
Danforth et al.ll201ll ) rather than the standard CalCOS flat field. 



Surprisingly, we find that the severity of the non-Poissonian noise in coadded COS 
data depends on the coaddition algorithm employed. In particular, we have found that the 
process of interpolating individual COS exposures onto a new wavelength scale (a neces- 
sary step in coadding multiple exposures with only partially overlapping wavelength ranges) 
can introduce additional non-Poissonian noise to the coaddition process. We have tested 
several algorithms, including linear, quadratic, cubic spline, and nearest neighbor interpo- 
lation schemes, and found that data coadded using nearest neighbor interpolation has the 
least severe non-Poissonian noise. The coadded data points in Figure (U for which we found 
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rj(x 



oc x 



0.37 



see Table [T]), utilize nearest neighbor interpolation. We find 77 oc x - 33 for 



data coadded using linear interpolation^] and rj oc x 36 for data coadded using quadaratic 
or cubic spline interpolation. We also find that the "xldsum" files generated by CalCOS 
have rj oc x 0,33 , identical to the behavior we find for linear interpolation. 

Regardless of its origin, non-Poissonian noise has an appreciable effect on the limiting 
equivalent width achievable with COS data. In the following subsections we calculate equiv- 
alent width limits for unresolved (Section 12. ip and spectroscopically resolved (Section 12.21) 
absorption lines that acc ount for the non-Poisson ian n oise in the d ata, and compare them 
to the limits derived by iGhavamian et al.l ( 120091 ) and iKrisd (120 111 ), which assume purely 
Poissonian noise. 



2.1. Unresolved Lines 

If we substitute the results of Equation [3] into Equation [Q, we find 

which reduces to Equation [2] in the Poissonian limit. The function w(x) has a unique 
minimum for all values of < j3 < 1 (see Eqn. H]) if f c is monotonic with x, as is the case for 
any reasonable LSF. This implies that a single value, x op t, exists that sets the most stringent 
limit on WW Physically, x opt represents the largest number of pixels that can be summed 
over a line profile befor e the noise contributions begin to outweigh the flux contributions 



( IGhavamian et al.l 120091 ) . 



Figure H] shows how w(x) varies with x for coadded data (solid line) and individual 
exposures (dot-dashed line), as compared to Poissonian data (dashed line) at a wavelength 
of 1250 A. Non-Poissonian noise causes x opt to be smaller and the minimum value of w(x) 
to be larger for coadded data and individual exposures than in the Poissonian case. Both 
of these effects are caused by the measured (S/N) x being smaller than Poisson statistics 
predict. 

Table [2] illustrates how non-Poissonian noise affects the limiting equivalent width of 
unresolved absorption lines in COS G130M and G160M spectra. We list the FWHM, x op t, 



1 Our coaddition routine originally used a linear interpolation scheme, but has been modified to use 
nearest neighbor interpolation instead as a result of this testing. Consequently, we recommend that 
users who downloaded previous versions of our coaddition code upgrade to the latest version (v2.0) at 
http : / / casa. Colorado . edu/$\sim$danf orth/costools .html 
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and limiting 3cr equivalent width of the iGhavamian et al.l (120091 ) and iKrisa (1201 if ) LSFs as 
a function of wavelength for the G130M and G160M gratings, as well as the corresponding 
values for the best-fit Gaussian LSF from pre-launch thermal vacuum testing. For each 
LSF, we tabulate th ese quantities assuming purely Poissoni an noise for c omparison with the 
values in Table 1 of IGhavamian et al.l (120091 ) and Table 3 of iKrissI ( 120111 ). which we are able 
to repr oduce. Ther e is a negligible difference in FWHM and x opt between the two LSFs 
but the IKrissI (120111) LSF produces equiv alent width limits that are ~ 2% larger than those 
derived from the IGhavamian et al.l (120091 ) LSF. 



For the IKrissI (120111 ) LSF we also tabulate the FWHM, x opt , and limiting 3cx equivalent 
width values for individual COS exposures and coadded COS data. These values allow us 
to quantify the trends that we illustrated in Figure [2J for unresolved lines, x opt is ~ 15% 
smaller and Wii m (3cr) is ~ 15% larger for individual exposures than Poissonian data, and 
x op t is ~ 25% smaller and Wn m (3a) is ~ 20% larger for coadded data than Poissonian data. 
These trends hold for the Gaussian LSF as well but their magnitude is slightly smaller due 
to the absence of the extended wings of the measured COS on-orbit LSF. Finally, we note 
that x opt for individual exposures is always ~ 0.2 pixels larger than the FWHM of the LSF, 
and x opt for coadded data is always ~ 0.6 pixels smaller than the FWHM of the LSF (i.e., 
x opt ~ 7.4 pixels for individual exposures and x opt — 6.5 pixels for coadded data; see Table|2]). 



2.2. Spectroscopically Resolved Lines 

The equivalent width formalism of Equation [5] can be extended to spectroscopically 
resolved absorption lines by modifying the fractional area of the line profile contained within 
the region of integration, f c (x). While for unresolved lines the line profile chosen was simply 
the instrumental LSF, the line profile of resolved features is the convolution of the LSF with 
the intrinsic shape of the absorption line. Generically, the intrinsic absorption line profile is 
described by a Voigt function that may have multiple velocity components; however, since 
we are only considering our ability to detect very weak features when calculating Wn m it is 
appropriate to assume a single-component Gaussian line profile. 

Convolving the instrumental LSF with the intrinsic line shape introduces a generic 
dependence on the Doppler parameter, b, to f c (i.e., f c = f c (x,b)). One consequence of this 
dependence is that calculating Wn m requires assuming a 6-value for the absorption line of 
interest. Determining an appropriate b- value is most straightforward if one is setting a limit 
on a line of the same ionic species as one of the lines detected in the spectrum (e.g., if one 
is trying to set a limit on the weaker line of a doublet when t he st ronger line is detected). 



Furthermore, as discussed in Section 121 and IGhavamian et al.l (120091 ) . the wings of the COS 
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LSF are wavelength dependent, introducing a wavelength dependence to f c for COS data 
such that f c = f c (x, A, b). 

Table[3]lists the limiting equivalent width of resolved absorption features in COS G130M 
and G160M spectra assuming a Gau ssian li ne pr ofile with Doppler parameter, b. This Gaus- 



sian profile was convolved with the iKrissI ( 120111 ) LSF and we have tabulated the Gaussian 



and convolved FWHMs, as well as x opt , f c {x pt) and Wii m (3cr) for Poissonian data, individ- 
ual exposures, and coadded data. For sufficiently large 6-values (>20 kms -1 ) the convolved 
FWHM is always ~ 2-3 pixels larger than the Gaussian FWHM; for smaller b- values the line 
is only marginally resolved and the convolution has a larger effect. Thus, the convolution 
acts to increase the apparent b- value of a line by approximately 1.5-2 pixels (~ 4 kms -1 at 
1450 A; see Eq. [7] and [S]) if the effects of this convolution are not accounted for. 

The effect of non-Poissonian noise on x opt for resolved lines is nearly identical to its 
effect on x opt for unresolved lines (see Section |2~T| Table |2]), but its magnitude decreases by 
~ 3% between b = 10 kms -1 and b = 200 kms -1 as the shape of the convolved line profile 
becomes more and more Gaussian. The effect on Wii m (3cx) is striking, however, with the 
limits for individual exposures increasing from 120% to 150% of the Poissonian values, and 
the limits for coadded data increasing from 130% to 190% of the Poissonian values, as the 
Doppler parameter increases from b = 10 kms -1 to b = 200 kms -1 . 



3. Empirical Scaling Relations 

We are now interested in determining the significance level of an absorption feature 
detected by the COS medium resolution gratings. We can do so by manipulating Equations [T] 
and replacing the limiting equivalent widths (Wn m ) defined therein with an arbitrary 
equivalent width, W\: 

N a {x, A, b) = (S/N)^ — (6) 

= Pitt fc{x,X,b) 

AX x 



AX w(x,X,b)' 

where we have now explicitly listed the dependence of f c on wavelength and Doppler param- 
eter and w is defined as in Equation [5J While Equation E] is algebraically simple, f c depends 
on the details of the LSF, making the equation cumbersome to use in practice. In particu- 
lar, for sufficiently large W\, our assumption of an intrinsically Gaussian line shape with a 
single velocity component (see Section 12. 2j) breaks down and there is no way of determining 
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a unique value of f c . However, it has been our experience that observers are most interested 
in determining whether an absorption feature is detected above a modest given significance 
threshold (3a, 4cr, etc.); for lines that are clearly detected above this threshold (i.e., where 
our assumption about the intrinsic line shape typically breaks down), accurately determining 
the exact value of the significance level of the detection is generally not important. 

To maximize the significance level of an observed line, on e should integr a te oy er a 
discrete region x opt pixels wide and evaluate w at x = x opt . In iDanforth et al.l ( 1201 ll ) we 
made the simplifying assumption that the absorption profile's deviations from Gaussianity 
were small, such that x opt equals the Gaussian FWHM and fAxopt) ~ 0.761, to reduce 
Equation [6] to a more tractable form (see Eq. 4 of IDanforth et al.l 1201 ll ). For coadded data, 
this assumption is justified if b > 30 kms -1 , where x opt and f c (x opt ) both differ from their 
assumed values by < 10% (see Table E]). 

Here we search for parameterizations of x opt and f c {x op t) that are valid for a broader 
range of 6-values, allowing us to reduce Equation [6] to a scaling relation between easily 
me asured quant ities while properly accounting for both the on-orbit COS LSF as defined 
by iKrissI (120111 ) and the non-Poissonian noise characteristics of the data as described in 
Section[2j To do so, we calculate a; opt and / c (x op t) at 13 wavelengths and 15 b- values spanning 
the ranges A = 1150-1750 A and b = 1-1000 kms" 1 , respectively. We parametrize both x opt 
and f c {x opt ) as functions of the Doppler line width, Xb, or the number of pixels required to 
span a Doppler parameter, b, at wavelength A: 



x b 



b 

c AA' 



(7) 



where c is the speed of light and AA is the dispersion of the spectrum. This width and the 
Gaussian FWHM listed in Table [3] are straightforwardly related by 



Gaussian FWHM = 2Vln2x fe « 1.665 x b . 



In IDanforth et al.l (120111 ) , we assumed that the optimal integration width equalled the 
Gaussian FWHM (Eq. [8]), but it's clear from Table [3] that the relationship is more compli- 
cated. We find that x opt is well fit by the function 



Xopt \Xb) 



a x x b + b x x b 01 , 



(9) 



where the best-fit values of the coefficients a\-Ci for Poissonian data, individual COS expo- 
sures, and coadded COS data are listed in Table [H The top panels o f Figures [3H5] show our 
fits to x opt (xb) (solid line) and the Gaussian approximation used by Danforth et al.l ( 2011 ) 
(dashed line) plotted over the values derived from the LSF model of iKrisd (120111 ). The 
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Gaussian approximation is nearly identical to our fits for x\, > 40 pixels (6 > 100 kms -1 ), 
but severely underestimates x opt for Xb < 10 pixels (b < 25 kms -1 ). For very small values of 
the Doppler width (x b < 1 pixel), Equation overestimates the value of x opt ; for these small 
values of Xb, we recommend using x opt = 6.5 pixels for coadded data, x opt = 7.4 pixels for 
individual exposures, and x opt = 8.7 pixels for Poissonian data, as shown by the dot-dashed 
horizontal lines in Figures [3H5J 

The bottom panels of Figures [3H5] show that / c (x op t) asymptotic ally approaches values 



that are slightly different than the constant value that we assumed in iDanforth et al.l (j201lh 



(dashed line). We parameterize f° pt (xb) = f c (x opt , X,b) with the function 

f^\x b )=a 2 -b 2 e-^ c \ (10) 

which is shown by the solid lines in the bottom panels of Figures [3H5J For Xb > 10 pix- 
els, G130M data (shown with asterisks in Fig. [3]) and G160M data (diamonds) agree very 
well, but there is a larger amount of scatter at smaller values of x b , which correspond to 
b < 20 kms -1 where absorption lines are becoming increasingly unresolved. Despite this 
scatter, Equation [10] does a reasonable job of modelling f° pt {x b ) for all values of Xb stud- 
ied, particularly when compared to assuming a constant value. The best-fit values of the 
coefficients a 2 -C2 for Poissonian data, individual COS exposures, and coadded COS data are 
listed in Table [TJ 

Substituting Equations 191 and [TU1 into Equation [6] and defining N° pt (xb) = N a (x opt , A, b), 
we find 

Wx f opt 

(x b = S/No pt ^^- (11) 

Wi 1 

= (S/N) 1 ^ 



AA w(x, 



opt J 



where (S/N) opt is the signal-to-noise measured at a resolution of A opt = x opt x AA. We 
have tested the accuracy of Equation [11] over the entire COS G130M + G160M bandpass 
(1150-1750 A) at 6-values of 1-1000 kms -1 and find t hat it predic ts significance levels that 



deviate from the values derived from the LSF model of lKrisd (120111 ) by < 10% over the entire 
parameter space, and < 5% when b > 10 kms -1 . 

While it is convenient to parameterize x opt and /° pt , and therefore w(x opt ), as a function 
of the single parameter x b , it is still useful to examine how wavelength and Doppler parameter 
affect these quantities individually. Figure [6] shows how f° pt and w(x opt ) vary as a function 
of wavelength for coadded data evaluated at Doppler parameters of b = 10, 30, 50, 100, 
and 200 kms -1 . Both quantities have a mild wavelength dependence but vary substantially 
with Doppler parameter. For modest values of the Doppler parameter (10 kms -1 < b < 
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50 kms" 1 ), / c opt increases quickly as the b- value increases before saturating at f° pt ^0.75 
(see Fig. [3] and Tabled]). Interestingly, w(x opt ) shows the opposite behavior — slow change 
for b < 50 kins -1 and a rapid increase for higher values. We only present the explicit 
wavelength dependence of /° pt and w(x opt ) for coadded data because we expect it will be 
the most applicable for observers, but the analogous behavior of individual exposures and 
Poissonian data can be straightforwardly derived from Equations (3HU and I7T-TT01 using the 
appropriate values from Table [H 

To facilitate the proper calculation of significance levels and equivalent width limits for 
coadded COS G130M and G160M data, we have written two short IDL functions that evalu- 
ate Equations ITl-fTTl The first, COS_SlGLEVEL, returns the significance level of an absorption 
feature in coadded data given an observed equivalent width, observed wavelength, estimated 
Doppler parameter, and either (S/N) opt or (S/N) 1 as inputs. The second, cos_ewlim, re- 
turns a limiting equivalent width for coadded data given a significance level, observed wave- 
length, estimated Doppler parameter, and either (S/N) opt or (S/N) 1 as inputs. These func- 
tions are available on the COS Tools website at|http : //casa. coloradoTedu/$\sim$danf orth/costools A 



3.1. The Benefits of Coaddition 

In the previous Sections we have shown that the presence of non-Poissonian noise has 
an appreciable effect on the significance of absorption features detected in COS FUV data 
and that the non-Poissonian noise is more severe for coadded data than individual exposures. 
These results may call into question the benefits of coadding COS data in the first place. 
We address these questions here. 

One obvious benefit to obtaining data at multiple FP-POS positions and/or central 
wavelength settings and c oadding to obtain th e final data product is that detector artifacts, 



such as gain sag regions (jSahnow et al.ll201ll ). that are not properly removed by standard 



CalCOS processing, will be present at slightly offset wavelengths in the individual expo- 
sures. These features move around in wavelength space precisely because they are fixed in 
detector (i.e., pixel) space, allowing for straightforward identification and removal, if neces- 
sary. 

To assess the more subtle tradeoffs of coadding data, we compare the noise properties 
of coadded data and individual exposures by rearranging Equation [IT] to find 

w ( x< o& dd ) (S/N)f add A^; ndiv yy™ add AA indiv 

M:Cp d t iv ) ~ (S/N)f div iV™ add W» ndiv AA coadd ' ^ 
where w^xl^ 11 ) and w(x™ dlv ) are w(x opt ) evaluated for coadded data and individual expo- 



- 12 - 



sures, respectively, and the other terms are defined analogously. This ratio is plotted as 
a function of x^ in Figure which shows that w(x™f d ) is approximately 6% larger than 
w ( x opt lv ) f° r unresolved lines, and their ratio increases as x^- 053 (dashed line) for resolved 
lines (xjj > 10 pixels). 

Equation [12] can be interpreted in several useful ways. If one assumes that the signal- 
to-noise, significance level, and dispersion of coadded data and individual exposures are the 
same, then Equation [T2] quantifies how much larger limiting equivalent widths for coadded 
data are than limiting equivalent widths for individual exposures. Similarly, if one assumes 
that the signal-to-noise, equivalent width, and dispersion of coadded data and individual 
exposures are the same, then Equation [12] quantifies how much larger the significance level 
for individual exposures are than the significance level for coadded data. Finally, and most 
germane to this discussion, if one assumes that the dispersion of coadded data and individual 
exposures are the same, then Equation [12] quantifies the improvement in signal-to-noise ratio 
required for an absorption feature of a given equivalent width to have the same significance 
level in coadded data and individual exposures. 

If two exposures of equal exposure time are coadded, then one would (naively) expect 
that (S/N)™ add = V^(S/N)i ndiv , in which case the increased signal-to-noise achieved through 
coaddition compensates for the increased non-Poissonian noise introduced by the coaddition 
process when xj, < 700 pixels (b < 1500 kms -1 ). We have made no efforts to test how 
signal-to-noise increases as a function of exposure time for COS data, but we note that even 
if coadding two equal-length exposures only increases the signal-to-noise per pixel by a factor 
of 2 37 w 1.29 (chosen to match the exponent /3 found when characterizing rj(x) for coadded 
data; see Tabled]), then coadding will still increase the overall sensitivity for Xb < 125 pixels 
(b < 300 kms" 1 ). 

Therefore, coadding two equal-length exposures will improve the limiting equivalent 
width for all but the broadest absorption lines. The limiting equivalent widths of even the 
broadest lines can be improved by adding more than two exposures together since we do 
not find any dependence of r)(x) on the number of exposures contributing to the coaddition 
(Section|2j). Since most COS FUV observations consist of four or more exposures per grating, 
coaddition will almost always improve the sensitivity of a dataset to weak absorption features. 

This happy result is largely a consequence of the individual exposures themselves being 
non-Poissonian. If we compare coadded COS data to Poissonian data we find that w(xl'^ dd ) 
is approximately 20% larger than w(x^l sson ) for unresolved lines, and their ratio increases as 
2.0.142 £ Qr reso i vec i lines. Thus, if individual COS exposures were Poissonian then coadding 
two equal-length exposures would only improve the limiting equivalent width for absorption 
lines with b < 25 kms -1 , assuming (S/N)^ oadd oc t 1//2 . 



13 



Conclusions 



Ghavamian et al.l ( 120091 ) and iKrisa (1201 if ) showed that the extended wings of the on- 
orbit LSFs of the COS medium resolution gratings lead to limiting equivalent widths for 
unresolved absorption lines that are ~ 30-40% higher than those derived from a Gaussian 
LSF, assuming Poissonian noise properties. We have demonstrated that COS G130M and 
G160M data have non-Poissonian noise characteristics and that the smoothed signal-to-noise 
ratio, (S/N)^, is proportional to x 0A2 for individual exposures and x ' 37 for optimally coadded 
data (see discussion in Section [2]), rather than x 1//2 as Poissonian statistics predict. When this 
non-Poissonian noise is accounted for, the limiting equivalent width of unresolve d absorption 
lines increases by ~ 15-20% over the iGhavamian et al.l ( 120091 ) and iKrissI ( 120111 ) values. 



We have also extended the equivalent width formalism of IGhavamian et al.l (120091 ) to 
spectroscopically resolved lines and find that the non-Poissonian noise properties of COS data 
have an even larger effect on resolved lines than unresolved lines. For a Doppler parameter 
of b = 10 kms -1 , the limiting equivalent width for coadded data is 30% larger than the 
corresponding Poissonian value; this percentage grows to 90% larger than the Poissonian 
value for b = 200 kms -1 . This effect is caused by the smoothed signal-to-noise differing 
more and more from the Poissonian value as the integration width increases, and emphasizes 
the importance of properly accounting for the non-Poissonian noise properties of COS G130M 
and G160M data. 

To that end, we have derived empirical scaling relations (Eqs. EHTT]) that allow easy 
calculation of the optimal significance of a detected absorption feature, N° pt , as a function 
of the Doppler line wi dth, X h - Th ese relations have been compared to the values derived 
from the LSF model of IKrissI (120111 ) in the wavelength range A = 1150-1750 A and Doppler 
parameter range b = 1-1000 kms" 1 and found to differ by < 10% over the whole parameter 
range, and < 5% when b > 10 kms -1 . IDL functions that use our parameterizations to 
determine the significance level of an observed feature or set an equivalent width limit on a 
nondetection have been developed and made publicly available (see Section [3]). 

While the non-Poissonian noise in coadded data is more severe than in individual ex- 
posures, the fact that the individual exposures are themselves non-Poissonian means that 
carefully coadding data can lead to significant reductions of the limiting equivalent width 
(i.e., improvements in sensitivity). Taking the non-Poissonian noise properties of both indi- 
vidual exposures and coadded data into account, we find that the signal-to- noise per pixel 
of coadded data need only increase by 6% over the signal-to-noise in individual exposures 
to reduce the limiting equivalent width of unresolved lines. For resolved lines, the necessary 
improvement in signal-to- noise scales approximately as x° 053 for Xb > 10 pixels (see Fig. [7]). 
Thus, a signal-to-noise improvement of < 25% is required for coadded data to have higher 
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sensitivity than individual exposures for lines with b < 150 kms . 

This work was supported by NASA grants NNX08AC146 and NAS5-98043 to the Univer- 
sity of Colorado at Boulder for the HST/COS project. We also thank T. Tripp and J. Tum- 
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Table 1. Best-Fit Coefficients 



Poissonian Individual Coadded 
Parameter Data Exposures Data Eqn 



r] 0.08 ±0.01 0.15 ±0.02 (4) 

(5 1/2 0.42 ±0.01 0.37 ±0.01 (4) 

ai 1.954 ±0.003 1.748 ± 0.002 1.605 ±0.002 (9) 

6i 6.9 ±0.1 5.8 ±0.1 5.1 ±0.1 (9) 

ci 0.17 ±0.01 0.22 ±0.01 0.25 ±0.01 (9) 

a 2 0.837 ±0.001 0.783 ± 0.001 0.743 ±0.001 (10) 

b 2 0.183 ±0.005 0.185 ±0.005 0.185 ±0.004 (10) 

c 2 10.5 ±0.4 11.0 ±0.4 11.6 ±0.4 (10) 



Table 2. Limiting Equivalent Widths for Unresolved Absorption Lines 







Ghavamian LSF (Poisson) 


Kriss LSF (Poisson) 


Kriss LSF (Indiv. Exp.) 


Kriss LSF (Coadded) 


Grating 


A 


FWHM 




W Um (3a) 


FWHM 


•^opt 


Wu m (3(7) 


■^opt 


Wum(3<r) 


^opt 


W Um (3a) 




(A) 


(pixels) 


(pixels) 


(mA) 


(pixels) 


(pixels) 


(mA) 


(pixels) 


(mA) 


(pixels) 


(mA) 


G130M 


Gauss 


6.5 


7.6 


9.9 


6.5 


7.6 


9.9 


6.6 


11.1 


5.8 


11.7 


G130M 


1150 


7.5 


9.2 


13.7 


7.4 


9.2 


14.1 


7.6 


16.2 


6.8 


17.2 


G130M 


1200 


7.3 


9.0 


13.5 


7.3 


9.0 


13.9 


7.6 


15.9 


6.6 


17.0 


G130M 


1250 


7.3 


9.0 


13.4 


7.2 


9.0 


13.8 


7.4 


15.8 


6.6 


16.8 


G130M 


1300 


7.2 


8.8 


13.3 


7.1 


8.8 


13.6 


7.4 


15.6 


6.6 


16.6 


G130M 


1350 


7 


8 fi 


1 3 1 


7 n 


8 6 


1 3 5 


7 9 


1 5 4 


fi 4 


16 3 


G130M 


1400 


7.0 


8.4 


13.0 


6.9 


8.4 


13.3 


7.2 


15.2 


6.2 


16.1 


G160M 


Gauss 


6.5 


7.5 


12.1 


6.5 


7.5 


12.1 


6.5 


13.7 


5.9 


14.4 


G160M 


1450 


7.3 


9.0 


16.1 


7.4 


9.0 


16.5 


7.5 


18.9 


6.7 


20.1 


G160M 


1500 


7.1 


8.8 


15.9 


7.2 


8.7 


16.3 


7.4 


18.6 


6.5 


19.8 


G160M 


1550 


7.1 


8.7 


15.9 


7.1 


8.7 


16.2 


7.4 


18.5 


6.5 


19.7 


G160M 


1600 


7.1 


8.7 


15.8 


7.0 


8.7 


16.1 


7.4 


18.3 


6.5 


19.5 


G160M 


1650 


7.0 


8.5 


15.6 


7.0 


8.5 


15.9 


7.2 


18.1 


6.4 


19.3 


G160M 


1700 


6.9 


8.5 


15.4 


7.0 


8.3 


15.7 


7.2 


17.9 


6.4 


19.0 


G160M 


1750 


7.0 


8.5 


15.3 


7.0 


8.5 


15.6 


7.2 


17.7 


6.4 


18.8 



Note. — All equivalent width limits assume S/N = 10 per pixel. The two line spread functions are detailed in 
Krissl (j201lh . 
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Table 3. Limiting Equivalent Widths for Resolved Absorption Features 



Poisson Indiv. Exp. Coadded 

Gaussian Convolved 



b 


FWHM 


FWHM 




/cOopt) 


Wi im (3<7) 


^opt 


fc{x op t) 


Wi im (3a) 


^opt 


/cOopt) 


Wim(3(T) 


(km s ) 


(pixels) 


(pixels) 


(pixels) 




(mAJ 


(pixels) 




(mAJ 


(pixels) 




(mAJ 












G130M (A 


= 1250 A) 












10 


7.0 


10.9 


13.2 


0.706 


15.5 


11.4 


0.649 


18.4 


10.0 


0.605 


20.1 


20 


13.9 


17.5 


21.1 


0.758 


18.1 


18.1 


0.696 


22.5 


16.3 


0.652 


25.2 


30 


20.9 


24.3 


29.1 


0.784 


20.6 


25.3 


0.726 


26.3 


22.9 


0.683 


30.1 


40 


27.9 


31.2 


37.3 


0.801 


22.8 


32.5 


0.743 


29.8 


29.5 


0.699 


34.6 


■50 


34.8 


38.1 


45.3 


0.809 


24.9 


39.9 


0.755 


33.1 


36.1 


0.710 


38.8 


75 


52.2 


55.3 


66.0 


0.824 


29.5 


58.0 


0.768 


40.5 


53.0 


0.726 


48.5 


100 


69.6 


72.6 


86.5 


0.830 


33.5 


76.2 


0.775 


47.0 


69.6 


0.733 


57.3 


150 


104.5 


107.1 


127.4 


0.835 


40.4 


112.5 


0.781 


58.6 


103.1 


0.740 


73.0 


900 


1 39 3 


141.6 


1 68 3 


837 

U.UU 1 


46.4 


149.1 


784 


68 Q 


1 36 4 


749 


87 












G160M (A 


= 1550 A) 












10 


7.0 


10.7 


12.9 


0.715 


18.5 


11.1 


0.659 


21.9 


10.0 


0.616 


23.9 


20 


14.1 


17.2 


20.6 


0.758 


22.0 


17.8 


0.700 


27.2 


16.2 


0.660 


30.5 


30 


21.1 


24.1 


29.0 


0.786 


25.1 


25.2 


0.729 


32.1 


22.7 


0.684 


36.7 


40 


28.2 


31.2 


37.3 


0.801 


28.0 


32.5 


0.744 


36.5 


29.6 


0.702 


42.4 


50 


35.2 


38.1 


45.6 


0.811 


30.6 


39.9 


0.754 


40.6 


36.3 


0.712 


47.7 


75 


52.8 


55.7 


66.4 


0.824 


36.3 


58.4 


0.768 


49.8 


53.3 


0.727 


59.8 


100 


70.4 


73.2 


87.2 


0.830 


41.3 


76.9 


0.776 


58.0 


70.2 


0.733 


70.6 


150 


105.6 


108.1 


128.5 


0.835 


49.8 


113.7 


0.781 


72.3 


104.0 


0.740 


90.1 


200 


140.8 


143.0 


170.1 


0.837 


57.1 


150.6 


0.784 


85.0 


137.9 


0.743 


107.5 



Note. — All equivalent width limits assume S/N = 10 per pixel. All convolutions use the COS line spread function detailed in lKriss 



(|201lh . 
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Fig. 1. — The ratio of the smoothed signal-to- noise to the signal-to-noise per pixel, r)(x), as 
a function of smoothing length, x. For purely Poissonian noise, r](x) = x 1 ^ 2 (dashed line). 
We find rj(x) ~ x°' 42 + 0.08 (dot-dashed line) for individual exposures and rj(x) ~ x°' 37 + 0.15 
(solid line) for coadded data (see Tabled]), indicating the presence of non-Poissonian noise. 
Only smoothing lengths larger than the FWHM of the COS on-orbit LSF (dotted vertical 
line; see Table [2]) contribute to the fits. 
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Fig. 2. — The function w(x) (see Eqn. [5]) evaluated at A = 1250 A for coadded data (solid 
line) and individual exposures (dot-dashed line), compared to the Poissonian value (dashed 
line). Non- Poissonian noise leads to a smaller value of x opt and a larger minimum for the 
function. The dotted vertical line shows the approximate FWHM of the COS on-orbit LSF, 
which is approximately equal to x opt for individual exposures (see Table |2]). 
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Fig. 3. — The optimal integration width, x opt (top panel), and the fractional area of an 
absorption line contained within the optimal integration width, f° pt (bottom panel), as a 
function of Doppler line width, Xb, evaluated for coadded COS FUV data. Simulated G130M 
data are shown with asterisks and G160M data with diamonds. The solid lines show our 
best fits to the data (Eq. [9] and [10] Table [TJ, which hold over the entire range of Xb studied 
although the data for /° pt show lar ger scatter when Xb < 10 pixels. The dashed lines show 
the Gaussian relations assumed in iDanforth et al.l (120111 ) . The dotted vertical line shows 
the approximate FWHM of the COS on-orbit LSF (see Table |2]), which roughly demarcates 
unresolved and resolved absorption lines. The dot-dashed horizontal line in the top panel 
shows our recommended value of x opt = 6.5 pixels for a; 6 < 1 pixel (see text). 
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Fig. 4. — Same as Figure El but evaluated for individual COS FUV exposures. The dot- 
dashed horizontal line in the top panel shows our recommended value of x opt = 7.4 pixels 
for xi, < 1 pixel (see text). 



-23 - 




0.1 1.0 10.0 100.0 1000.0 

x b [pix] 

Fig. 5. — Same as Figure El but evaluated for purely Poissonian data. The dot-dashed 

horizontal line in the top panel shows our recommended value of x opt = 8.7 pixels for 
Xb < I pixel (see text). 
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Fig. 6. — The dependence of /° pt (top panel) and w(x opt ) (bottom panel) on wavelength 
for coadded data evaluated at Doppler parameters (from top to bottom in both panels) of 
b = 200, 100, 50, 30, and 10 kms -1 . G130M data are shown with asterisks and G160M data 
with diamonds. The values predicted by our empirical scaling relations (Equations [9] and 
TUl) are indicated by a dashed line over the G130M wavelength range and a dot-dashed line 
over the G160M wavelength range. 
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Fig. 7. — The ratio w(x^ dd ) /w(x^ d ") (see Eq. [12]) . plotted both as a function of Doppler 
line width, Xb (bottom axis), and as an approximate function of Doppler parameter, b (top 
axis, for which the labelled values are only strictly valid when A = 1450 A). The dotted 
vertical line shows the approximate FWHM of the COS on-orbit LSF (see Table [2]), which 
roughly demarcates unresolved and resolved absorption lines. The dashed line shows the 
function x° 053 , which is an excellent fit to the plotted ratio for resolved lines. 



